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Abstract 

We present results on rapidity and transverse momentum distri- 
butions of inclusive charm quark production in heavy-ion collisions at 
RHIC, including the next-to-leading order, O(a^), radiative correc- 
tions and the nuclear shadowing effect. We find the effective, nuclear 
K-factor to be K(y) ~ 1.4 for | y |< 3 in the rapidity distribution, 
while 1 < K(pt) < 3 for lGeV < pt < 6GeV in the pt distribution. 
We incorporate multiple parton scatterings in our calculation of the 
fraction of all central events that contain at least one charm quark 
pair. We obtain the effective A-dependence of the charm cross sec- 
tions. Finally, we comment on the possibility of detecting the quark- 
gluon plasma signal as an enhanced charm production in heavy-ion 
collisions at RHIC. 



The main goal of the future heavy-ion colliders, such as RHIC, is to study 
the properties of nuclear matter under extreme conditions, and in particular 
to search for the formation of a new state of matter, the quark-gluon plasma 
[1]. With the assumption that high energy nuclear collisions lead to the 
thermalized system, RHIC energies are sufficiently large to produce very 
dense matter, of the order of several GeV / fm 3 , well above the critical density 
necessary to create the quark-gluon plasma [2]. The problem of finding a 
clean, detectable signal for this phase is currently one of the most challenging 
theoretical problems. In the last few years, the possible signals, such as 
thermal photons, dileptons and J/ty suppression, have been investigated in 
detail and found to have yields comparable in magnitude with those expected 
from a simple extrapolation of hadronic collisions [3] . Recently, open charm 
production in heavy-ion collisions has been proposed as an elegant method 
for probing the possible formation of the quark-gluon plasma [4]. However, 
in order to determine whether the enchanced charm production at RHIC 
can be unambiguously interpreted dbS db SI gnal of QGP, one needs to have 
control of the other sources of charm production. The standard competing 
process is charm production through the hard collisions of partons inside the 
nuclei. In addition to being relevant as a background for the signal of quark- 
gluon plasma, this type of charm production is useful tool for studying the 
perturbative aspects of strong interactions and for determining the nuclear 
screening/shadowing effect on the gluon distribution in a nucleus. 

In this letter we present results of calculations of the rapidity and trans- 
verse momentum distributions of inclusive charm quark production in Au-Au 
collisions at RHIC, including the 0(af) radiative corrections and the nuclear 
shadowing effect. We determine the size of the gluon contribution to the 
charm production at RHIC. We also calculate the fraction of central and in- 
elastic events at RHIC that will contain charm quarks which satisfy unitarity 
constraints by properly taking into account multiple nucleon scatterings. In 
our calculations we use the latest set of two-loop evolved parton densities in 
a nucleon obtained from global fits of data from deep inelastic lepton-nucleon 
collisions [5] . We take the nuclear shadowing effect in the quark distribution 
from the recent measurements in deep inelastic lepton- nucleus collisions [6]. 
We assume that the amount of shadowing present in the gluon distribution 
is the same as in the quark. To illustrate the importance of next-to-leading 
order contributions and the nuclear shadowing effects, we determine the effec- 
tive (i.e. nuclear) K-factor defined as the ratio of the particular distribution 
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to the leading-order distribution without any nuclear effects. We show that 
this K-factor is very different than the one in hadronic collisions, and that 
in general, it can not be approximated by a constant. We also determine 
the effective A-dependence of the charm cross sections. Finally, we comment 
on our results for the number of charm quarks produced in central rapid- 
ity region in central Au-Au collisions in the context of quark-gluon plasma 
signatures at RHIC. 

In perturbative QCD, the inclusive cross section for charm production 
in nuclear collisions is obtained by convolution of parton densities in nuclei 
with a hard scattering parton cross section [7]. In our calculation, for the 
parton cross sections, we include leading-order subprocesses, 0(a 2 ), such as 
q + q — > Q + Q and g + g — > Q + Q, and next-to-leading order contributions, 
0(a 3 s ), such &sq + q-*Q + Q + g,g + q— >Q + Q + g,g + q-*Q + Q + q and 
g + g^Q + Q + g. The double differential inclusive distribution of charm 
production in central AA collisions can be written as 

dN c ^ , , da c , . 

= T AA (0)— (1) 



cPpTdy d 2 p T dy ' 

where the double differential inclusive cross section is given by 

J^- = 1" / dM(Xa , Qi )F ?( Xi , Q 2 ) -^gff - * , (2) 

and Taa(O) is the nuclear overlapping density at zero impact parameter, 
F t A (x, Q 2 ) is the parton structure function in a nucleus, x a and Xb are the 
fractional momenta of the incoming partons, s is the parton-parton cm. 
energy (s = x a Xbs). The parton differential cross section calculated to 0(af) 
can be written as [8] 

4^r = - h ^ + 7^-Ah ( 3 ) 

d?p T dy s 1,3 2ns 2 J 

Expressions for the functions h\°j and h^j can be found in Ref. 8. The 
coupling constant a s (Q 2 ) that appears in Eq. (3) is given by 



a s (Q 2 ) = ^ oA l - 6(153 -™ lngVA2 ], (4) 

(33 - 2N f ) In % L (33-2A^ / ) 2 lng 2 /A 2 
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where Q 2 is the renormalization scale, A is the QCD scale parameter and 
Nf is the number of flavors. We take the factorization scale in the structure 
functions to be 2m c and we consider the renormalization scales Q = m c and 
Q = 2m c . We do not evolve the structure functions below Q = 2m c , because 
for Q 2 < 8.5GeV 2 their behavior is not very well known [9]. For the mass of 
the charm quark we use m c = 1.5GeV. 

The double differential distribution for inclusive charm quark production 
in inelastic Au-Au collisions can be obtained from Eq. (1) by replacing 
T AA (0) with A 2 /af r £ lastic , where af n A elastic « 4irR 2 A . 

To obtain the number of nucleon-nucleon collisions per unit of transverse 
area at fixed impact parameter, we consider the nuclear overlapping function 
[10] 

T AA {b) = J d 2 hT A {\ h \)T A (\ b- h |), (5) 

where the nuclear profile function, T A (b), is the nuclear density integrated 
over the longitudinal size, i.e. 

T A (b) = [°° dzp A (VWT^) (6) 

For nuclear density we use the Woods-Saxon distribution [11] given by 

The density and the nuclear overlapping function are normalized so that 
/ d 3 rp(r) = A and / d 2 bT AA (b) = A 2 . For central collisions the overlap- 
ping function can be approximated by T AA (0) = A 2 /irR A 2 , which gives 
T Au - Au (0) = 30.7771&- 1 . 

The nuclear parton distribution, if nucleons were independent, would be 
given as A times the parton structure function in a nucleon. However, at 
high energies, the parton densities become so large that the sea quarks and 
gluons overlap spatially and the nucleus can not be viewed as a collection 
of uncorrelated nucleons. This happens when the longitudinal size of the 
parton, in the infinite momentum frame of the nucleus, becomes larger than 
the size of the nucleon. Partons from different nucleons start to interact 
and through annihilation effectively reduce the parton density in a nucleus. 
When partons inside the nucleus completely overlap, there reach a saturation 
point. Motivated by this simple parton picture of the nuclear shadowing 
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effect and taking into account the A 1 ^ 3 dependence obtained by considering 
the modified, nonlinear modifying factor to the Altarelli-Parisi equations with 
gluon recombination included, the parton structure function in a nucleus can 
be written as [12] 



R(x,A) 



F z A {x,Q 2 ) 
AF % N (x,Q 2 ) 



= < 
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(8) 



where Ff is the parton structure function in a nucleon, x n = l/(2r p m p ), xa 
is a saturation point (xa = l/(2i?^m p )), m p is the mass of the proton, r p 
is the radius of a proton and Ra is the radius of the nucleus. The only free 
parameter is a constant D which can be determined by fitting the data. 

Comparison of the shadowing function with all deep inelastic lepton- 
nucleus data on the ratio F^(x,Q 2 )/F^(x,Q 2 ) [6] indicate that Eq. (8) 
has much steeper x-dependence, especially for 0.002 < x < 0.1, the region of 
relevance to charm production at RHIC and LHC energies. In addition, the 
onset of saturation for the ratio of structure functions for Xe to Deuterium 
is observed at values of x about order of magnitude smaller than those pre- 
dicted by the parton recombination model [6]. Thus, it is not surprising that 
even the best fit of Eq. (8) to the data overestimates the observed shadowing 
effect by about 15%. Consequently charm production calculated with this 
shadowing function would be underestimated by about 40%. 

In our calculation we use the shadowing function that has recently been 
proposed as the best fit of EMC, NMC and E665 data [11] and is given by 
[13] 

( a-i — Q>aX Xn < x < 0.6 

(9) 



R(x,A) 



[a 3 - a 4 x ) 1+kqAai{1/x _ 1/xo) x<x 



The values for the parameters k q , ai, a2, «3, «4 and x$ can be found in 
Ref. 13. It is interesting to note that for x ~ Xq = 0.14, this shadowing 
function has a form similar to the shadowing function motivated by the 
parton recombination picture. The main difference is the onset of saturation, 
which occurs at smaller value of x than xa of Eq. (8) and much more 
gradually. In addition, the A-dependence is much weaker, namely A 0A . 

By integrating Eq. (2) and multipling by Ty^(0) (for central collisions) 
or by A 2 /<j{^ (for inelastic collisions), we obtain rapidity distribution, trans- 
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verse momentum distribution and the total cross section for the inclusive 
charm quark production in AA collisions. We compare our results for the 
total charm cross section with the low-energy hadronic and nuclear data [14], 
for the energy range 20GeV < y/s < hhGeV . We find very good agreement 
with all the data for the choice of scale Q = m c . When Q = 2m c is used, 
our cross sections are slightly smaller than the measured values. Detailed 
comparison of our results with low-energy measurements of differential dis- 
tributions and the total cross section for charm production will be published 
separately in Ref. 15. 

Here we present our results on rapidity and transverse momentum distri- 
butions for the inclusive charm production in Au-Au collisions at RHIC. In 
our calculation we use the two-loop evolved parton structure functions, MRS- 
SO [5] with A 5 = lAOMeV. By using two other sets of structure functions, 
MRS-DO and MRS-D- [5], we find that theoretical uncertainty due to the 
choice of the nucleon structure function is only 10%. This is not surprising 
because the average x-value probed with charm production at ^/i = 200GeV 
is about 10 -2 , which is still within the range of x for which there is deep in- 
elastic lepton-nucleon scattering data. Another theoretical uncertainty is the 
choice of the renormalization scale. We find about 40% lower cross sections 
when we use the scale Q = 2m c instead of Q = m c . Since the low-energy data 
seem to be better described with the choice of scale Q = m c , we will present 
most of our results with that particular choice . Our results for differential 
distributions (da/dy, da/dpx, da/dxp and da / dpTdy\ y=0 ) and for the total 
charm cross section in heavy-ion collisions at LHC, (including the theoretical 
uncertainties), will be included in Ref. 15. 

Our result for the rapidity distribution of inclusive charm production in 
central Au-Au collisions at RHIC is presented in Fig. la (solid line). We 
also show the rapidity distribution when nuclear shadowing is not included 
(dotted line) and the leading-order results with shadowing (short-dashed line) 
and without shadowing (long-dashed line). We note the the shape of the 
rapidity distribution does not seem to be affected by the next-to-leading 
corrections or by the nuclear shadowing effect. In Fig. lb) we show our 
results for two different choices of scale, illustrating the uncertainty due to 
this particular choice. We find that, in the central rapidity region, the number 
of charm quark pairs produced per unit rapidity in central Au-Au collisions 
at RHIC is 0.6 for the scale Q = 2m c and 0.9 for the scale Q = m c . 

In hadronic collisions, one usually defines the "K-factor" as a measure of 



5 



the size of higher-order corrections. Here we define the effective K-factor for 
nuclear collisions as a ratio of the particular distribution to the leading-order 
distribution without any nuclear effects. In Fig. 2 we present our results 
for the K-factor. We show that in hadronic collisions K-factor is about 2 
(squares), while the nuclear K-factor is about 1.4 (circles) in the central 
rapidity region < 3). This is due to the fact that the nuclear shadowing 
effect effectively suppresses production of charm quarks by about 30%. 

In Fig. 3 we present our results for the transverse momentum distribution 
of the charm quark produced in Au-Au collisions at RHIC (solid line). We 
find that both higher-order correlations and the nuclear effects change the 
shape of pt distribution. The nuclear shadowing effect is much stronger at 
low pt (about 40% effect), while at px = 6GeV it reduces the cross section 
by only few percent. The next-to-leading order corrections give a factor of 
1.7 increase at low pt and about factor of 3 at pt = QGeV. These two effects 
together result in effective K-factor increasing from 1 at px — lGeV to 3 at 
Pt = 6GeV. At large pt, where nuclear shadowing effects are negligible, we 
expect K-factor to approach its hadronic value. We find similar behavior of 
the K-factor for xf distribution, namely its strong dependence on xf [15]. 

By integrating differential distributions over the phase space we obtain 
the total number of charm quark pairs produced. For central (inelastic) 
Au-Au collisions we get about 4 (1) charm quark pairs produced per event. 

To obtain the effective A-dependence of the total inclusive charm cross 
section in nuclear collisions, defined as <J^ t arm = we use the total 

charm cross section in hadronic collisions at y^s = 200GeV, o v c p = 0.18m6. 
We find that (3 = 1.27 for central collisions and (3 = 1.94 for inelastic colli- 
sions. 

To be able to determine the fraction of central or inelastic events which 
contain at least one charm quark pair, we need to consider the semiclassical 
probability of having at least one parton-parton collision at fixed impact 
parameter, 1 — e~ TAA ^ a ° , where e~ TAA ^ ac is the probability that there is 
no parton-parton scattering in Au-Au collision at impact parameter b. The 
fraction of events in Au-Au collisions that contain at least one charm quark 
pair is then given by 




/ <fb[l - exp(-T AA (b)a c )} 
Jd 2 b[l - exp(-T AA (b)af%)\ 



(10) 
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where Taa is given by Eq. (5) and and a c is the integrated charm cross 
section [] 

l l partons 

a c = dx a dx b [Fi/AiPa, Q 2 )F j/A (x b , Q 2 )a id (s, m 2 c , Q 2 ) (11) 

s x a s i,j 

The parton cross section aij(s,m 2 ,Q 2 ) has been calculated to the order 
O(a^) and can be written as [17] 

^M^ 2 ) = ^/,( P |), (12) 

where 

f id (p, Q 2 /m 2 c ) = f®(p) + Ana s (Q 2 )[fQ(p) + f$(p) In {Q 2 /m% (13) 

The functions //°- ,and f^j, are given in Ref. 17. 

To determine the fraction of all central events that contain at least one 
charm quark pair we integrate Eq. (9) over the small range of impact pa- 
rameter, i.e. < b < 0.1 fm. We find that about 98% of central events will 
contain at least one charm quark pair. For inelastic collisions we integrate 
Eq. (9) over all impact parameter and find this fraction to be 38%. Note 
that the integrated charm cross section in Eq. (9) includes multiple inde- 
pendent parton-parton scatterings which means multiple charm quark pair 
production. 

To conclude, we have presented the complete next-to-leading order cal- 
culation of the differential and total inclusive cross sections for the charm 
production in Au-Au collisions. We have shown that at RHIC energies, 
both higher-order contributions and the nuclear shadowing effect are large 
and can not be neglected. In the central region of the rapidity distribu- 
tion, the higher-order contributions increase the cross section by a factor of 
2, while the nuclear shadowing effect result in additional decrease of about 
30%. These two effects together result in nuclear K-factor of about 1.4. On 
the other hand, in case of the px distribution, the K-factor changes from 1 at 

1 The total cross section for heavy-quark production in hadronic collisions has been 
previously calculated using the Eq. (11) and found to be in agreement with UA1 data on 
bottom production at y/s — 630GeV [16]. 
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Pt = lGeV to 3 at pr = QGeV. The effective A-dependence for the charm 
production in the central (inelastic) collisions is found to be A 1 ' 27 (A 194 ). 
By properly incorporating multiple parton scatterings we have shown that 
about 98% (38%) of all central (inelastic) events at RHIC will contain charm 
quark pairs. We have found that at RHIC energies the dominant subprocess 
for charm production is gg fusion (about 95%). Therefore, future measure- 
ments of charm production in p-p, p-A and A-A collisions at RHIC energies, 
in addition to being a test of perturbative QCD, could provide valuable in- 
formation about presently unknown A-dependence and x-dependence of the 
nuclear shadowing in the gluon density. 

We would like to emphasize that our calculation of charm production 
includes next-to-leading order corrections and is performed in the standard 
parton model with the use of factorization theorems for hard scatterings in 
pQCD. As a consequence of the truncation of the perturbative series, our 
results are sensitive to the choice of the renormalization and factorization 
scale. We find this uncertainly to be about 40%, which is roughly an esti- 
mate for the size of the higher-order corrections. We evolve the factorization 
scale down to Q — 2m c . The behavior of the structure functions at scales 
below 2m c are not very well known [9]. Furthermore, we do not perform our 
calculation with renormalization scale below m c , because the perturbative 
calculation becomes unreliable. The nuclear shadowing effect is incorporated 
in a simple way, namely by using the parametrization of the recent data on 
shadowing of the quark densities [6] with the assumption that the shadowing 
effect for gluons is the same. The shadowing of the gluon density, obtained 
perturbatively by solving the modified Altarelli-Parisi equation for the struc- 
ture functions [13] is less than the observed effect for quarks [6]. Further the- 
oretical work along these lines is necessary to have a better understanding 
of parton densities in the nucleus. Another approach would be to consider 
the role of quantum-mechanical interference in the heavy-ion collisions [18]. 
The most interesting result of this novel approach is that both the shadow- 
ing and the antishadowing effects are obtained without any modification of 
the structure functions. We expect future high precision measurements of 
deep-inelastic p-A scatterings and the Drell-Yan production to be able to 
test both the standard parton model picture and the quantum-mechanical 
interference effect providing an important information about the origin of 
the nuclear shadowing effect. 

Finally, we make a remark on the possibility of detecting a signal for the 
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formation of quark-gluon plasma via enchanced charm production in heavy- 
ion collisions at RHIC. We have found that 0.6 — .9 open charm quark pairs 
per unit rapidity (in the central region) will be produced in central Au- 
Au collisions via hard parton-parton scatterings. Even though our results 
seem to indicate that unrealistically high initial temperature of the quark- 
gluon plasma [4] is needed to overcome this size of the background, further 
theoretical and experimental work is needed in quantitative understanding of 
the nuclear shadowing effect, especially of the gluon density, before definite 
conclusion can be made. 
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Figure Captions 



Fig. 1. a) Rapidity distribution of inclusive charm quark production in 
Au-Au collisions at RHIC, calculated to next-to-leading order (LO+NLO) 
including nuclear shadowing (NS) (solid line), without NS (dotted line), 
only leading-order (LO) without NS (long-dashed line), and with NS (short- 
dashed line) , b) The same calculation as above but for two different choices 
of the scale, Q 2 = ml (solid line) and Q 2 = Ami (dotted line). 

Fig. 2. The effective K-factor, K = (^)/(^)lo, for the distributions 
presented in Fig. 1. The K-factor for the next-to-leading order (LO+NLO) 
distribution without NS (squares), including NS (circles) and for the leading- 
order (LO) distribution with NS (diamonds). Note that K w 2 for hadronic 
collisions (squares). 

Fig. 3. The same as in Fig. la) but for the transverse momentum 
distribution of inclusive charm quark production. Curves are labeled as in 
Fig. la). 

Fig. 4. The effective K-factor for the transverse momentum distributions 
presented in Fig. 3. Labeling is the same as in Fig. 2. 
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